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ABSTRACT 


Under  suitable  programs  of  mechanical  or  thermal  loading,  many  solid  mate¬ 
rials  are  capable  of  undergoing  phase  transformations  from  one  crystal  structure 
to  another.  The  austenite-martensite  transformation  that  occurs  in  a  variety  of 
metallic  alloys,  including  the  so-called  shapeTmemory  materials,  provides  an  exam- 
pie.  The  present  paper  represents  an  effort  to  model  coupled  thermo-mechanical 
effects  in  the  macroscopic  response  of  solids  that  arise  from  the  occurrence  of  phase 
transformations.  A  Helmholtz  free  energy  potential  is  constructed  to  describe  the 
thermo-fnechanical  response  of  the  hypothetical  material  to  be  considered  here.  As 
a  function  of  shear  strain,  the  potential  is  non-convex  in  a  certain  range  of  tem¬ 
perature;  this  feature  is  essential  for  the  modeling  of  phase  transformations.  Apart 
from  some  general  preliminary  considerations  pertaining  to  finite  thermoelasticity, 
the  analysis  is  carried  out  in  the  context  of  a  simple  problem,  idealized  from  an  ex¬ 
periment,  in  which  an  annular  cylinder  is  deformed  to  a  state  of  radically  symmetric, 
finite  anti-plane  shear  in  the  presence  of  differing  inner  and  outer  surface  tempera¬ 
tures.  After  constructing  all  radially  symmetric  weak  solutions  involving  at  most  a 
single  surface  of  discontinuity  of  strain  or  temperature  gradient,  we  determine  the 
implications  for  quasi-static  motions  of  the  second  law  of  thermodynamics.  The 
therm<£hiechanical  phase  transformation-induced  hysteresis,  residual  deformation 
and  stress  relaxation  effects  exhibited  by  this  model  are  discussed.  The  results  con¬ 
cerning  creep  rate  as  predicted  by  the  present  model  are  in  qualitative  agreement 
with  the  laboratory  observations.  Finally,  the  shape-memory  effect  as  predicted  by 
the  present  thermo-mechanical  model  in  the  setting  of  finite  anti-plane  shear  has 
been  simulated  numerically;  a  number  of  pictures  selected  from  the  simulation  are 
included. 
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1.  Introduction 

Many  solid  materials  are  capable  of  undergoing  phase  transformations  from 
one  crystal  structure  to  another  under  suitable  programs  of  mechanical  or  ther¬ 
mal  loading.  The  austenite-martensite  transformation  that  occurs  in  a  variety  of 
metallic  alleys,  including  the  so-called  shape-memory  materials,  provides  a  well 
known  example.  Experimental  observations  of  the  behavior  of  such  materials  are 
discussed,  for  example,  in  [1].  The  development  of  models  for  the  macroscopic 
thermo-mechanical  response  of  such  materials  represents  an  active,  current  area 
of  study  in  the  continuum  mechanics  of  solids.  Examples  of  such  investigations 
may  be  found  in  [2]- [10]. 

The  problem  of  modeling  the  macroscopic  response  of  a  solid  undergoing  a 
phase  transformation  is  an  inherently  nonlinear  one.  Most  of  the  studies  reported 
in  [2]- [10]  were  carried  out  within  the  framework  of  the  purely  mechanical  theory 
of  finite  elastostatics;  they  are,  therefore,  appropriate  only  for  the  modeling  of 
load-induced  phase  transformations  taking  place  at  a  given  temperature.  In  these 
studies,  the  interface  between  two  material  phases  emerges  as  a  surface  of  strain 
discontinuity  in  the  relevant  equilibrium  field.  Not  all  elastic  potentials  -  or 
strain  energy  functions  -  are  capable  of  sustaining  equilibrium  fields  exhibiting 
such  singular  surfaces.  Those  that  do  have  this  capability  correspond  to  what 
have  been  called  non-elliptic  elastic  materials,  because  they  lead  to  displacement 
equations  of  equilibrium  that,  although  elliptic  for  infinitesimal  deformations,  lose 
their  ellipticity  for  a  certain  range  of  finite  strain. 

When  surfaces  of  strain  discontinuity  occur  in  equilibrium  elastic  fields,  quasi¬ 
static  motions  -  regarded  as  one-parameter  families  of  equilibrium  states  -  become 
dissipative.  Dissipation  occurs  in  the  sense  that  the  rate  /ork  of  the  external 
forces  acting  on  a  portion  of  the  body  no  longer  coincides  w  1  Ihe  rate  of  increase 
of  stored  energy.  Their  difference  may  be  represented  as  an  integral  over  the  slowly 
moving  surface  bearing  the  strain  discontinuity  of  a  fictitious  “driving  traction”  - 
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or  “Eshelby  force”  [11]  -  multiplied  by  the  velocity  of  the  moving  singular  surface. 
In  this  isothermal,  purely  mechanical  setting,  the  disparity  between  the  rate  of 
external  work  and  the  rate  of  storage  of  energy  may  thus  be  regarded  as  arising 
because  work  must  be  done  to  move  the  phase  boundary. 

Equilibrium  boundary  value  problems  in  the  theory  of  non-elliptic  elastic 
materials  suffer  a  massive  breakdown  in  the  uniqueness  of  solution.  Quasi-static 
motions  in  such  materials  may  visit  any  of  the  infinitely  many  available  equi¬ 
librium  states  consistent  with  the  boundary  conditions.  In  order  to  render  the 
macroscopic  response  of  such  materials  determinate,  it  is  necessary  to  supplement 
the  constitutive  description  of  the  material  by  a  relation  that  specifies  the  “kinet¬ 
ics”  of  the  phase  transformation,  or  the  rate  at  which  particles  transform  from 
one  phase  to  another.  The  incorporation  of  such  a  kinetic  relation  into  the  theory 
has  been  shown  to  predict  isothermal  macroscopic  response  that  is  in  qualitative 
agreement  with  experiments,  at  least  in  the  one-dimensional  theory  of  tension 
experiments  [13,1]. 

The  present  paper  represents  an  effort  to  model  coupled  thermo-mechanical 
effects  on  the  macroscopic  response  of  solids  that  arise  from  the  occurrence  of 
phase  transformations.  Apart  from  some  general  preliminary  considerations  per¬ 
taining  to  finite  thermoelasticity  in  Chapter  2,  our  analysis  is  carried  out  in  the 
context  of  finite  anti-plane  shear,  which  is  the  simplest  setting  in  which  to  study 
finite  deformations.  Chapter  3  is  devoted  to  the  construction  of  a  Helmholtz  free 
energy  potent  ial  intended  to  describe  the  thermo- mechanical  response  of  the  hypo¬ 
thetical  material  to  be  considered  here.  This  material  is  designed  to  model  phase 
transitions  in  somewhat  the  same  sense  as  the  van  der  Waals  model  describes  the 
gas-liquid  transition.  As  a  function  of  strain,  the  potential  is  non-convex  in  a 
certain  range  of  temperature. 

Motivated  by  the  possibility  of  large  strains  arising  from  phase  transforma¬ 
tions  in  the  earth’s  mantle,  Sammis  and  Dein  [14]  carried  out  a  simple  laboratory 
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experiment  designed  to  show  that  creep  could  indeed  be  generated  by  the  occur¬ 
rence  of  a  phase  transformation  induced  by  both  thermal  and  mechanical  loads. 
Their  experiment,  which  is  described  in  detail  in  Chapter  4,  is  one  that  lends  it¬ 
self  to  modeling  by  finite  anti-plane  shear,  and  Chapters  4  through  8  are  devoted 
to  the  formulation  and  analysis  of  such  a  model.  In  particular,  the  implications 
for  quasi-static  motions  of  the  second  law  of  thermodynamics  are  determined  in 
Chapter  5,  and  construction  of  a  model  kinetic  relation  is  undertaken  in  Chapter 
6.  Chapter  7  deals  with  the  phase  transformation-induced  hysteresis,  residual 
deformation  and  stress  relaxation  effects  exhibited  by  the  model.  The  results 
presented  in  Chapter  8  concerning  creep  rate  as  predicted  by  the  present  model 
of  the  Sammis-Dein  experiment  are  in  qualitative  agreement  with  the  laboratory 
observation  reported  in  [14]. 

Chapter  9  illustrates  the  shape-memory  effect  as  predicted  by  the  present 
simple  thermo-mechanical  model  in  the  setting  of  finite  anti-plane  shear. 
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2.  Preliminaries  on  continuum  mechanics 

2.1  Global  and  local  versions  of  balance  laws 

Consider  a  body  B  that,  in  a  reference  configuration,  occupies  a  region  R 
in  three-dimensional  Euclidean  space.  A  motion  of  the  body  on  a  time  interval 
L  =  [to, ti]  is  characterized  by  a  one-parameter  family  of  invertible  mappings 
y(-,t)  :  R  i — ►  Rt,  with 

y  as  y(x,  t)  =  x  +  u(x,  t)  for  (x,t)  6/2x1.  (2.1) 

We  assume  that  the  deformation  y,  or  equivalently  the  displacement  u,  is  con¬ 
tinuous  with  piecewise  continuous  first  and  second  derivatives  on  /txL.  Let 
F(x, f)  =  Grady(x,t)  stand  for  the  deformation  gradient  tensor,  and  assume 
that  for  each  t  E  L,  and  x  E  R,  the  Jaccobian  determinant  of  mapping  (2.1)  is 
positive: 


J(x,t)  =  de*F(x,t)  >  0.  (2.2) 

Let  p(x)  denote  the  mass  density  of  B  at  the  point  x  in  the  reference  region 
R,  b(x,  t)  the  body  force  per  unit  mass,  and  <r(x,t)  the  nominal  stress  tensor 
during  the  motion.  At  each  t,  we  require  p(-)  and  b(-,<)  to  be  continuous  on  R, 
while  <r(-,  t)  is  to  be  piecewise  continuous  with  a  piecewise  continuous  gradient  on 
R.  Let  p*  (y,  t )  be  the  mass  density  of  B  at  the  point  y  in  the  region  Rt  occupied 
by  the  body  at  the  instant  t,  and  we  require  p*(-,  t)  to  be  piecewise  continuous  for 
each  t.  For  each  t  E  L,  the  balance  laws  for  mass,  linear  and  angular  momentum 
require 


anda  + 


(2.3) 


(2.4) 
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f  y  A  anda  +  f  y  A  pbdv  =  -j-  [  y  A  pixdv ;  (2.5) 

Jau  J d  jd 

respectively,  for  all  regular  subregions  D  C  R.  The  superposed  dot  denotes  the 
partial  derivative  with  respect  to  t  at  fixed  x. 

Next,  let  q(x,  t)  denote  the  nominal  heat  flux  vector,  £(x,i)  the  heat  supply 
per  unit  mass  and  e(x,  t)  the  internal  energy  per  unit  mass.  At  each  t,  we  suppose 
that  £(■,£)  is  continuous  on  R  and  that  q(-,<)  is  piecewise  continuous  with  a 
piecewise  continuous  gradient  on  R.  The  internal  energy  c(-,  •)  is  required  to  be 
piecewise  continuous  with  piecewise  continuous  first  derivatives  on  R  x  L.  The 
first  law  of  thermodynamics  requires  that  at  each  instant  t, 

I  <m  •  u da+  I  pb  •  udv  +  /  q  •  n da  +  /  p£dv 
JdD  JD  J8D  JD  ( o  p\ 

d  r  t2'6^ 

=  jt  J^p(e  +  u-u/2)dvt 

for  every  regular  subregion  D  C  R.  Finally,  let  0(x,  £)  denote  the  absolute  tem¬ 
perature  and  rj(x,t)  the  entropy  per  unit  mass.  At  each  t,  we  assume  that  6(-,  t) 
is  continuous  with  a  piecewise  continuous  gradient  on  R,  while  is  assumed 
to  be  piecewise  continuous  with  piecewise  continuous  first  derivatives  on  R  x  L. 
The  rate  of  entropy  production  in  a  region  D  in  12  is  defined  by 

r (£;  D)  =  ^  [  prjdv  -  [  q  •  n/0da  -  [  pt/ddv.  (2.7) 

Jd  JdD  Jd 

The  Clausius-Duhem  version  of  the  second  law  of  thermodynamics  requires 


r(£;  D)>  0  D  C  R,  teL. 


At  a  fixed  instant  t,  localization  of  the  balance  laws  (2.3)  -  (2.6)  and  the 
inequality  (2.8)  at  a  point  x  at  which  all  fields  are  smooth  yields  the  following: 
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p*J  -  p , 

Div  a+pb  =  pii, 

<tFt  =  F<rT,  (2.9) 

<r-F  +  Div  q  +  p£  =  pc, 

Dxu  (q/0)  +  p£/0  <  prj. 

On  the  other  hand,  suppose  that  S(t)  is  a  regular  surface  in  R  at  time  t  across 
which  some  or  all  of  the  thermo-mechanical  quantities  suffer  jump  discontinu¬ 
ities.  Localization  of  (2.3)  -  (2.8)  at  a  point  x  on  S(t)  yields  the  following  jump 
conditions: 


\p*A  =  o, 

[<m]  +  p[i iR  =  0, 

1  (210) 

•  u]  -f  [p(e  +  -u  •  u)]Vn  +  [q  •  nj  =  0, 

W^n  +  lq  •  n/^1  <  0, 

where 

V„=V.n,  (2.11) 

and  V  =  V(x,  t)  is  the  velocity  of  the  moving  surface  S(t)  at  the  point  x.  The 
unit  normal  n  on  the  singular  surface  S(t)  is  chosen  such  that  Vn  >  0;  if  Vn  >  0, 
the  positive  side  of  S(t)  is  the  side  into  which  V  (and  therefore  n)  points.  If 
g(x,t)  denotes  a  generic  field  quantity  that  jumps  across  S(t),  we  write  [p(x,  t)J  = 
p+(x,t)  —  g~(x,t),  where  p+(x,t)  and  g~(x,t)  stand  for  the  limiting  values  of  g 
at  the  point  x  on  S(t)  from  the  positive  and  negative  sides,  respectively. 

Conversely,  the  field  equations  (2.9)  together  with  the  jump  conditions  (2.10) 
imply  the  global  balance  laws  (2.3)  -  (2.6)  and  (2.8). 

In  addition  to  the  jump  conditions  listed  in  (2.10),  one  can  derive  the  fol¬ 
lowing  kinematic  results  through  the  smoothness  requirements  imposed  on  the 
deformation  (2.1): 


[F]m  =  0,  [u]  =  -[F]V  on  S(t), 

where  m  is  any  vector  tangent  to  the  singular  surface  S(t). 


(2.12) 


2.2  Rate  of  entropy  production 

With  help  of  the  field  equations  (2.9)  and  the  jump  conditions  (2.10)  and 
(2.12),  Abeyaratne  and  Knowles  [12]  derived  the  following  alternate  representa¬ 
tion  for  the  rate  of  entropy  production  defined  in  (2.7): 


r(t;  D)  =  rloc(t;  D)  +  IWt;  D)  +  r,(t;  D),  (2.13) 

where 

rjoc(<;  D)  =  f  (pqQ  +  ff  t  -  pe)/0dv,  (2.14) 

Jd 

r con{t‘,D)=  I  q-GradO/O2  dv,  (2.15) 

Jd 

r*(t;  D)  =  f  J(<r++<r-).Fj  VjOda,  (2.16) 

Js(t)nD  * 

while  V*  is  the  Helmholtz  free  energy  per  unit  mass,  defined  by 

t)  =  e(x,  t)  -  0(x,  t)ij(x,  t)  on  R  x  L.  (2.17) 


In  (2.13),  the  total  rate  of  entropy  production  r(f;£>)  at  the  instant  t  in  the 
subregion  D  C  R  is  decomposed  into  three  parts:  arises  from  local  dissipation 

in  the  material  away  from  the  singular  surface;  is  the  entropy  production 
rate  due  to  heat  conduction;  finally,  r,  represents  the  entropy  production  rate 
due  to  the  moving  singular  surface  S(t).  No  constitutive  assumptions  are  made 
in  the  derivation  of  (2.13)  -  (2.16). 

The  driving  traction  f{x,t)  acting  on  the  surface  of  discontinuity  S(t)  is 
defined  by 
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/(x,  t)  =  P(x)[ip(x,  t)  1  -  - (<r+(x,  t)  +  <r  (x,  t ))  •  [F(x, t)], 

/or  x  G  S(t),  t  6  L. 


(2.18) 


Then,  from  (2.16)  and  (2.18),  r,(i;  D)  can  be  written  as 


r,(t;£>)=  f  fVjeda , 

JS(t)nD 


S(t)nD 

and  localization  of  (2.8)  at  point  x  on  S(t)  yields 


(2.19) 


/(x,t)Vrn(x,t)  >  0  for  x  €  S(i),  t  e  L. 


(2.20) 


2.3  Finite  thermoelasticity 

We  now  specify  that  the  body  at  hand  is  thermoelastic  and  assume  that  there 
are  a  free  energy  potential  and  a  heat-flux-response  function  q  such  that,  for 
every  thermodynamic  process, 


xJ>  =  *(F,0),  (2.21) 

q  =  q(F,0,  Grad  0),  (2.22) 

<t  =  p&f(F,0),  or  ffij  =  pd4?(T,d)/dFij,  (2.23) 

77  =  — ^g(F,  0),  or  tj  = —d4'(F,0)/d0.  (2.24) 

Objectivity  requires  that 

*(QF,  0)  =  F,  0 ),  q(QF,  0,  v)  =  q(F,  0 ,  v),  (2.25) 


for  every  nonsingular  tensor  F,  every  proper  orthogonal  tensor  Q,  every  positive 
scalar  0  and  every  vector  v. 
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Because  of  (2.17),  (2.21),  (2.23)  and  (2.24),  rjoc  defined  by  (2.14)  vanishes 
identically  for  thermoelastic  materials,  and  hence  localization  of  (2.8)  away  from 
the  singular  surface  S(t)  yields 


r  (•,  •,  v)  •  v  >  0,  for  any  vector  v.  (2.26) 

This  supplies  a  restriction  on  the  heat-flux-response  function  q. 

2.4  Isotropic  thermo-elastic  materials 

We  assume  from  here  on  that  the  body  is  thermo-mechanically  isotropic,  so 

that 


#(FQ,  S)  =  *(F,  6),  q(FQ,  0,  QTv)  =  QTq(F,  8,  v),  (2.27) 

for  every  nonsingular  tensor  F,  every  proper  orthogonal  tensor  Q,  every  positive 
scalar  0  and  every  vector  v.  Here  QT  denotes  the  transpose  of  Q.  It  follows  that 
the  free  energy  ^(F,0)  depends  on  F  only  through  the  deformation  invariants: 


*(F  ,0)  =  ^(/1,J2,/3,0), 

where  the  fundamental  scalar  invariants  h,h>h  are  defined  by 

h(C)  =  frC,  I2( C)  =  [(trC)2  -  tr(Ca)]/2, 

J3(  C)  =  detC,  C  =  FrF. 

In  the  presence  of  isotropy,  (2.23)  takes  the  form 


(2.28) 


(2.29) 


er  =  2pVIlF  +  2p# Ia(hl  -  FFT)F  +  2I3p^iaF-T,  (2.30) 

where  1  stands  for  the  identity  tensor  and  F“T  is  the  transpose  of  the  inverse  of 
the  nonsingular  tensor  F. 

We  assume  from  now  on  that  the  heat-flux-response  function  q  depends  on 
Grad 9  linearly,  i.e., 
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q(F,0,v)  =  K(F,0)v,  (2.31) 

for  every  nonsingular  tensor  F,  positive  scalar  9  and  vector  v.  Here  K  is  called 
the  heat  conductivity  tensor  and,  in  the  presence  of  isotropy,  K  takes  the  form 

K(F,  6)  =  fol  +  fa¥TF  +  02(FtF)2,  (2.32) 

where  fa  =  fa(h,I2tI3,9),  k  =  0,  1,  2. 

By  assuming  that  <f> i  and  fa  both  vanish  identically,  we  obtain  the  generalized 
Fourier  law: 


q(F,  9,  Grad  9)  =  fa(IuI2,  h ,  9)Grad9.  (2.33) 

fa  is  required  to  be  non-negative  by  (2.26). 

2.5  Thermo-mechanical  states  of  anti-plane  shear 

Suppose  .ft  is  a  cylindrical  region  and  let  X  =  (0;  ©i,  e2,  eg)  be  a  rectangular 
cartesian  frame  with  the  unit  vector  es  parallel  to  the  generators  of  the  cylin¬ 
der.  Consider  a  time-independent  thermo-mechanical  state  of  anti-plane  shear, 
in  which 


Vi=*i»  V2  =  X2,  y3  =  x3  +  u(xi,x2),  9  =  9(x  i,*2)»  (2.34) 

where  the  out-of-plane  displacement  u  and  the  temperature  9  are  both  defined  on 
the  cross-section  D  of  ft.  In  the  given  frame,  the  deformation  gradient  tensor  F 
has  components: 


Fa0  —  Fa3  =  0,  F3a  =  u,a  F33  =  1.  (2.35) 

Here  6ap  is  the  Kronecker  delta,  Latin  and  Greek  subscripts  have  the  respective 
ranges  1,  2,  3,  and  1,  2,  and  repeated  subscripts  are  summed  over  the  appropriate 
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range.  A  subscript  proceded  by  a  comma  indicates  partial  differentiation  with 
respect  to  the  corresponding  x-coordinate.  The  fundamental  invariants  7i,  J2,  J3 
are  given  by 


4  =  4  =  3-H72,  4  =  1,  (2.36) 

where  7  =  \Gradu\  =  y/u,a  u,a.  When  specialized  to  the  anti-plane  shear  (2.34), 
the  constitutive  statement  (2.30)  yields  the  components  of  nominal  stress  a  as 

crat 9  =  [2 &Ix  +  2(2  +  72)^/a  -  p]s<*i 9  “  2^/ati,a  u,p  , 

O-aS  =  [— 2^/j  +  p]u, a  ,  <r3a  =  2[^/x  +  §/a]u,a  ,  (2.37) 

^33  =  2#/,+ 4  §/3-p, 

where  we  have  written 


*  =  -b(h,h,e)  =  f>Mh,h,i,6),  (2.38) 


and  from  (2.36), 


71  =  72  =  3  +  72- 


(2.39) 


In  the  absence  of  body  forces,  the  momentum  balance  (2.9)2  yields 


[p  -  2 *h  -  2(2  -1-  72)^/al,a  +[2^/au,o  U  =  0»  a  =  1, 2,  (2.40) 

and 

[M  (7,  &)ui0  ],/3  =  0,  (2.41) 

where  the  arguments  7i  and  72  of  p  and  ^  in  (2.40)  are  given  by  (2.39),  and  M 
in  (2.41)  is  defined  by 
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M(7,0)  =  2  [#/,(/!, J2,0)  +  ^h(hJ2te)}Il=h=3+^.  (2.42) 

In  the  absence  of  heat  sources,  the  first  law  of  thermodynamics  (2.9)4  takes 
the  simple  form 

[*(%  »)«..], «  =  0,  (2.43) 

where  k  is  defined  by 

fc(7.0)  =  *>(3  +  72.3  +  72.l,«)-  (2.44) 

In  view  of  the  fact  that  p,  M  and  k  are  expressed  in  terms  of  u  and  9  through 
(2.38),  (2.42)  and  (2.44),  the  equations  (2.40),  (2.41)  and  (2.43)  actually  comprise 
four  differential  equations  for  two  unknowns  u  and  9,  and  hence  constitute  an 
over-determined  system  in  general.  Without  some  restriction  on  the  material 
properties,  one  would  expect  that  only  very  special  solutions  of  (2.41)  and  (2.43) 
would  also  satisfy  (2.40).  We  say  that  a  thermoelastic  material  is  capable  of 
sustaining  anti-plane  shear  if,  for  every  domain  D,  every  solution  u,  9  of  (2.41) 
and  (2.43)  also  satisfies  (2.40).  In  general,  an  isotropic  thermoelastic  material 
need  not  be  capable  of  sustaining  anti-plane  shear.  We  turn  next  to  a  sub-class 
of  such  materials  which  do  sustain  anti-plane  shear. 
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3.  A  special  class  of  thermoelastic  materials 

As  shown  by  Knowles  [15] ,  all  homogeneous,  isotropic  incompressible  hyper- 
elastic  materials  whose  elastic  potential  W  depends  only  on  the  first  fundamental 
scalar  invariant  Ji  are  capable  of  sustaining  anti-plane  shear.  Materials  in  this 
class  have  been  called  generalized  Neo-Hookean  materials  by  Gurtin  [10]  because 
the  Neo-Hookean  model  of  rubber  elasticity  is  included  among  them.  Generalized 
Neo-Hookean  materials  have  proved  to  be  useful  in  studying  nonlinear  effects  -  es¬ 
pecially  qualitative  ones  -  in  the  simplest  possible  setting:  that  of  finite  anti-plane 
shear;  see,  for  example,  [10],  [16],  [17]. 

Generalized  Neo-Hookean  materials,  being  incompressible,  do  not  permit  the 
useful  incorporation  of  thermal  effects  of  the  kind  essential  in  the  study  of  ther¬ 
moelasticity.  Jiang  and  Knowles  [18]  have  introduced  a  class  of  homogeneous, 
isotropic,  compressible  hyperelastic  materials  that  do  not  suffer  from  this  limita¬ 
tion  and  at  the  same  time  are  capable  of  sustaining  finite  anti-plane  shear.  The 
thermoelastic  counterpart  of  the  class  of  materials  studied  in  [18]  is  characterized 
by  the  Helmholtz  free  energy  per  unit  reference  volume  as  follows: 


*  (Ji.  /a,  h,0)  =  *(h ,  e)  +  *7l  (h ,  e)g(I3, 9)  +  g(I3 , 9 ),  (3.1) 

where  9  is  the  absolute  temperature,  and  4r,  f  and  g  are  functions  on  (0,  oo)  x 
(0,  oo),  with  ^  differentiable  three  times,  g  and  g  twice.  One  can  show  that  the 
necessary  and  sufficient  condition  for  a  material  of  type  (3.1)  to  be  capable  of 
sustaining  finite  anti-plane  shear  is  the  following 


<j(l,0)  =0,  p/s(l,0)  = -1,  gja(l, 9)  =  constant,  V  9  >  0.  (3.2) 

In  all  the  following  discussion,  we  choose  gi3(l,9)  to  be  zero,  which  implies  that 
the  undeformed  state  is  unstressed.  Without  loss  of  generality,  we  set  <7(1,0)  =  0 


and  thus  represents  the  Helmholtz  free  energy  density  associated  with  volume- 
preserved  deformations*. 


The  components  of  nominal  stress  in  (2.37)  reduce  to 

<Ta0  =  0,  <Ta3  =  <T3a  =  Affr,  0)u,a  ,  ^33  =  0,  OH  D,  (3.3) 

where  M  is  given  by  (2.42).  The  true  (or  Cauchy)  stress  r  is  given  in  terms  of  a 
and  F  by 


r  =  &FT  /J.  (3.4) 

From  (3.4),  (3.3)  and  (2.35),  one  can  determine  the  true  stresses  ry: 

rap  =  0,  Ta 3  =  r3a  =  M (7, 0)uta ,  r33  =  7 2M (7, 6) .  (3.5) 

We  introduce  the  shear  stress  response  function  f(7, 0)  by 

r(7 . 0)  =  (7, 5),  /or  (7, 0)  e  [0, 00)  x  (0, 00),  (3.6) 

then  the  magnitude  of  the  resultant  of  shear  stresses  ra3ea  is  then  given  by 

Vra3Ta3  =  r(7, 0).  (3.7) 

One  can  show  that  for  0  fixed,  the  equilibrium  equation  in  X3  .direction  (2.41) 
is  elliptic  [15]  for  any  out-of-plane  displacement  u  if  and  only  if 

Mh,0)dT(£e)  >  0  V7>0.  (3.8) 

or  equivalently, 

4,1(3  +  72.#)l'®;1(3+72,«)  +  27sf/,/.(3  +  72.«)]>0,  V7>0.  (3.9) 

*  Note  that,  $(/i,/2,/3,9)  =  if  h  =  1,  by  (3.2). 


Consider  a  deformation  of  a  body  composed  of  a  homogeneous,  isotropic 
thermoelastic  material  with  the  Helmholtz  free  energy  potential  ^(/i,/2,/3,0). 
If  Ai,  A2,  A3,  are  the  principal  stretches  of  the  deformation,  one  has 


h  =  A?  +  \\  +  Xl  I2  =  AjAj  +  AjA|  +  A3A*,  J3  =  Af  AjA*.  (3.10) 

Let 

*(Ai,A2,A 3,0)  =  (3.11) 

where  I\,  /2  and  h  are  expressed  in  terms  of  the  stretches  A»  by  (3.10).  It  follows 
from  (3.4),  (2.30),  (3.11)  and  (3.10)  that  the  principal  true  stresses  Ti,  r2  and  73 
are  given  by 


A  d^l 

= -^^(Ai.Aij.As,^),  {nosumoni).  (3.12) 

The  Baker-Ericksen  inequalities  require  that  the  greater  of  any  two  princi¬ 
pal  stresses  shall  correspond  to  the  greater  of  the  corresponding  two  principal 
stretches  when  the  stretches  are  distinct;  thus 


{ji  -  Tj ) (A*  -  A j)  >0  if  A<  /  Xjt  *,  j  =  1,  2,  3.  (3.13) 

For  the  material  described  by  (3.1),  one  finds  from  (3.1),  (3.10)  -  (3.12)  that  the 
Baker-Ericksen  inequalities  (3.13)  hold  [18]  if  and  only  if 


Vii(h,0)  +  ^fi1il{h,0)f{l3t6)  >  0,  for  0  <  I3  <  (Ji/3)3.  (314) 

When  specialized  to  the  anti-plane  shear  (2.34),  (3.14)  reduces  to 

^/i(3  +  T2  *0)  >  0,  for  7  >  0,  (3.15) 

i 


or  equivalently 
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M ( 7, 8)  >  0,  for  7  >  0.  (3.16) 

In  deformations  other  than  anti-plane  shear,  issues  such  as  the  form  taken  by 
the  ellipticity  conditions  for  the  material  characterized  by  (3.1)  are  more  difficult 
to  analyze.  The  implication  of  idealizations  in  (3.1)  is  thus  not  fully  clear.  One 
idealization  that  is  apparent  from  (3.2)  is  the  fact  that  the  reference  configuration 
is  unstressed  at  any  temperature. 

We  shall  be  concerned  with  a  special  case  of  the  materials  (3.1)  whose 
Helmholtz  free  energy  for  volume-preserving  deformations  is  given  by 

*(/l,  9)  =  -^(0)(/ l  -  3)[Ji  -  3  -  4/3(0)  v/TTr3  +  6a2)  +  fc(0),  (3.17) 

for  Ii  >  3,  8  >  0; 

with 

0(8)  =  AKa/A,)"^1-®/®0*3,  for  8  >  0;  (3.18) 

Here,  the  functions  y  and  h  axe  defined  and  twice  continuously  differentiable  on 
(0,  oo)  and  axe  required  to  satisfy 

y(8)  >  0,  h"(8)  <  0,  for  8  >  0;  (3.19) 

where  the  constant  do  is  a  charateristic  temperature  and  a,  0o,  m  axe  dimension¬ 
less  positive  constants  subject  to  the  restriction 

(0o/a)2  >  4/3,  (3.20) 

this  allows  the  material  to  achieve  some  unstressed  but  deformed  states  at  certain 
range  of  temperature  (see  eq.(3.28)  and  Fig.  3a)  and  this  feature  is  essential  to 
model  shape-memory  behavior  (see  Chapter  9). 
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For  convenience,  let 

*(7,  e)  =  tt(3  -(-  72, 0),  far  7  >  0,  0  >  0;  (3.21) 

by  (3.6),  (2.42)  and  (3.17),  the  shear  stress  response  function  f  at  constant  tem¬ 
perature  of  this  material  is  given  by 

f  (7,  0)  =  *,(7. «)  =  J»*(9)7(7a  -  30(9)7  +  3oJ)-  (3.22) 

In  this  special  case,  the  ellipticity  condition  (3.8)  reduces  to 

(72  -  30(0)7  +  3a2)(72  -  2 0(0)7  +  a2)  >  0,  V  7  >  0.  (3.23) 

From  (3.18),  one  can  see  that  (3.23)  fails  for  0  e  [0*,0*]  where 

9.  =  (1  -  i-)90,  9'  =  (1  +  1)0„.  (3.24) 

m  m 

Let  Ai,  A2  and  >13  be  the  sets  of  points  in  the  strain-temperature  plane 
defined  by 

Ax  =  {(7,  0)/  7  >  0  for  0  <  0«  orcd  0  <  7  <  7a/ (0)  /or  0  6  [0*,  0*]}, 

A2  =  {(7,*)/  7a/ (0)  <  7  <  7m (0)  /or  0  6  (0.,0*)},  (3.25) 

^3  =  {(7, 0)/  7  >  0  far  0  >  0*  aruf  7  >  7m(0)  /or  0  €  [0*,  0*]}, 
where  the  functions  of  temperature  7a/ (0)  and  7m(0)  are  defined  on  [0*,  0*]  by 

7tn(0)  =  m  +  VW)  -  “J.  7m(0)  =  0(9)  -  v/02(0)  -  “J.  (3.26) 

By  (3.21)  and  (3.17), 

{<0  on  A2\ 

=  0  on  dA2\ 

>0  on  A\  U  ^3. 


(3.27) 
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Here,  '£(•,0)  is  convex  for  0  G  (0,0*]  U  [0*,oo)  and  non-convex  for  0  G  (0*,0*). 
Correspondingly,  the  shear  stress  response  function  f(-,0)  is  monotonic  for 
0  €  (0,0.)  U  [0*, oo)  and  non-monotonic  for  0  G  (0*,0*)  (See  Figures  1  -  3). 
Because  these  features  of  the  free  energy  and  the  stress  response  change  at  the 
temperatures  0*  and  0*,  we  call  0*  and  0*  the  lower  and  upper  transition  temper¬ 
atures,  respectively. 

Noting  (3.17),  one  may  easily  show  that  the  special  version  of  the  Baker- 
Ericksen  inequalities  (3.15),  or  equivalently  (3.16),  fails  for  0  €  [0,,0*]  where 

S  L  H 4/3) 

*'  0  my  2ln((3o/a)  ’  (3.28) 

0*  =  20o  -  0*. 

Note  that  the  existence  of  0*  and  0*  is  guaranteed  by  (3.20).  By  (3.24),  (3.28) 
and  (3.20),  we  have 


0  <  0,  <  0.  <  0O  <  0*  <  0*.  (3.29) 

We  now  investigate  the  invertibility  of  the  shear  stress  response  function 
f(-,0).  By  (3.27),  there  exist  functions  rlf  T2  and  T3  such  that 

r  1  (t(7,  0) ,  0)  =7,  for  (7, 0)  G  Ai ; 

r2(r(7, 0),  0)  =7,  /or  (7, 0)  G  i2;  (3.30) 

r3 (r(7,  0) ,  0)  =7,  /or  (7, 0)  G  j43; 

with  continuities: 

rx(rM(9),9)  =  r2(r„(<?),0)  =  /or  9  6  [».,»•), 

r2 (rm(9),0)  =  r3(rm{«),»)  =  7m(«).  /or  0  e  [#..«*], 

(3.31) 

ri(f(7,0.),9.)  =  r,(f (7, ».),».),  for  7  >  a, 

r,(f  (7, «•).»*)  =  r,(f  (7,  *•),»•),  /or  a  >  7  >  0, 


Here  the  functions  rm(0)  and  rM{9)  are  defined  on  [0.,d‘]  by 
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Tm(0)  =  f(7m(0),  0),  tm(6)  =  f(yM(d),  0). 

From  (3.32),  (3.31)  and  (3.26),  one  can  see  that 

tm{&)  >  rm(0),  forde  (0„  0*); 
tm(8)  =  rm(0),  for  9  =  9,  or  0*. 

For  each  fixed  0  in  (0*,  0*),  there  is  a  number  re(9)  such  that 


(3.32) 


(3.33) 


rm(9)  <  rc(9 )  <  rM(9), 


(3.34) 


fV3(Tc(0),0) 

/  r(7, 9)  d7  =  Tc(0)[r3(rc(0),  9)  -  r^rc (9),  0)],  (3.35) 

J  Ti(tc{0),0) 

tc(0)  is  the  so-called  Maxwell  stress  at  the  temperature  9.  By  (3.34),  (3.33), 
(3.32),  (3.26),  (3.24)  and  (3.18), 


lim  rc(0)  -  rm(0*)  =  rM (0.)  =  f  (a,  0.), 
lim  rc(0)  =  rm(0 *)  =  tm{9*)  -  f(a,  9*). 


(3.36) 


To  characterize  a  thermoelastic  material  completely,  we  must  also  specify  the 
heat  conductivity  k  defined  by  (2.44).  From  here  on,  we  assume  that 


fc(7,0)  =  ki(9),  for  (7,0)  €  Au  i  =  1,  2,  3, 
where  fc»(-)  are  positive- valued  functions  and 


(3.37) 


*i(0)  *  fc 2(9 )  *  k3(9)  *  *i(0),  for  9  €  (0.,0*) 


(3.38) 
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4.  Equilibrium  of  an  annular  cylinder  in  anti-plane  shear 

In  [14],  Sammis  and  Dein  describe  a  simple  experiment  intended  to  illustrate 
the  role  of  phase  transformations  in  producing  large  strains  and  creep  (“superplas¬ 
ticity”)  in  solids.  Their  interest  in  this  issue  was  motivated  by  the  possibility  that 
superplastic  behavior  in  the  earth’s  mantle  might  arise  from  thermo-mechanically 
induced  phase  transformations. 

The  experiment  of  Sammis  and  Dein  involves  an  annular  cylindrical  crystal  of 
cesium  chloride  (see  Fig.4)  whose  outer  cylindrical  surface  is  clamped  and  main¬ 
tained  at  a  uniform  high  temperature.  The  inner  cylindrical  boundary  is  attached 
to  a  glass  tube  through  which  cold  water  flows  for  the  purpose  of  maintaining  the 
inner  wall  at  a  uniform  low  temperature.  A  weight  is  hung  from  the  end  of  the 
glass  tube  to  provide  an  axial  traction  on  the  inner  wall.  The  specimen  is  thus 
subject  to  both  a  mechanical  load  and  a  temperature  gradient.  For  outer  wall 
temperatures  below  a  certain  critical  value,  the  crystal  responds  in  a  conventional 
thermoelastic  manner,  with  the  inner  surface  displaced  axially  by  definite  amount 
that  depends  on  the  size  of  the  axial  load  and  the  temperature  gradient.  When 
the  temperature  at  the  outer  wall  exceeds  the  critical  value,  a  phase  transforma¬ 
tion  from  body-centered  to  face-centered  cubic  crystal  structure  is  initiated  at 
the  outer  wall,  and  the  interface  between  the  two  phases  moves  radially  inward. 
While  the  phase  boundary  is  moving,  the  weight  attached  to  the  central  glass 
cooling  tube  experiences  substantial  axial  creep,  a  phenomenon  that  would  not 
occur  if  the  response  were  that  of  a  conventional  thermoelastic  material.  Sammis 
and  Dein  measured  the  axial  creep  rate  at  the  inner  wall  as  a  function  of  time 
for  both  subcritical  and  supercritical  outer  wall  temperatures.  In  the  subcritical 
case,  very  low  creep  rates  were  observed,  consistent  with  expectations  based  on 
conventional  thermoelastic  response.  In  contrast,  much  larger  creep  rates  were 
observed  when  the  outer  wall  temperature  exceeded  the  critical  value. 


Our  next  objective  is  to  formulate  a  problem  designed  to  provide  a  simple 
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qualitative  model  for  the  experiment  of  Sammis  and  Dein[14].  We  model  the 
specimen  in  their  experiment  as  an  annular  cylinder  composed  of  a  thermoelastic 
material  in  the  class  of  materials  specified  in  the  preceding  chapter.  We  idealize 
the  experiment  itself  as  one  in  which  the  annular  cylinder  is  deformed  to  a  state  of 
radially  symmetric,  finite  anti-plane  shear  in  the  presence  of  differing  inner  and 
outer  surface  temperatures.  The  anti-plane  shear  is  produced  by  a  prescribed, 
uniformly  distributed  axial  traction  acting  on  the  inner  surface.  After  formulating 
the  corresponding  equilibrium  boundary  value  problem,  we  construct  all  radially 
symmetric  weak  solutions  involving  at  most  a  single  surface  of  discontinuity  of 
strain  or  temperature  gradient,  which  we  call  an  equilibrium  shock.  In  Chapters  5 
and  6,  we  discuss  quasi-static  anti-plane  shear  motions  of  the  annular  cylinder,  and 
in  Chapter  7  we  determine  the  corresponding  dissipative  macroscopic  response. 
A  qualitative  comparison  of  the  creep  rates  observed  by  Sammis  and  Dein  with 
the  results  predicted  by  the  present  model  is  carried  out  in  Chapter  8. 

Chapter  9  is  devoted  to  an  illustration  of  the  shape-memory  effect  in  the 
context  of  the  model  put  forward  here. 

4.1  Formulation  and  reduction  of  the  problem 

Let  the  open  reference  cross-section  D  be  the  region  bounded  by  two  concen¬ 
tric  circles  of  radii  a  and  b  (a  <  b;  see  Fig.  5).  The  temperatures  on  the  inner 
wall  and  outer  wall  of  the  cylinder  are  both  specified  constants.  Let  the  outer 
wall  be  fixed,  and  let  the  axial  shear  traction  acting  on  the  inner  wall  be  given. 
Thus,  the  boundary  conditions  are  given  by 

T 

9  —  0a,  oy 3  =  — - —  where  r  =  a,  (4.1) 

27ra 


0  =  6h, 


lt3  =  0 


where  r  =  b. 


(4.2) 
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Here  r,  <f>  and  23  form  a  cylindrical  coordinate  system,  and  the  constant  T  is  the 
resultant  force  per  unit  axial  length  acting  on  the  inner  wall.  We  consider  radially 
symmetric  anti-plane  shear  displacement  and  temperature  that  are  independent 
of  time: 


u  =  u(r),  0  —  0(r)  on  D  (4.3) 

We  suppose  the  body  to  be  composed  of  the  thermoelastic  material  described 
by  (3.1)  and  (3.17)  in  Chapter  3.  Hence,  in  the  absence  of  body  force,  the  in¬ 
plane  equilibrium  equations  (2.40)  are  satisfied  automatically  and  the  out-of-plane 
equilibrium  equation  (2.41)  may  be  integrated  to  give 

f(7(r),  0(r))  =  r  €  [a,  b],  (4.4) 

where  f  is  given  by  (3.6).  The  energy  equation  (2.43)  reduces  to 

M7(r)i0(r))““j“  —  2nr'  f  6  M’  <45> 

where  go  stands  for  the  constant  value  of  heat  flux  in  the  cylinder.  Using  (2.10) 
and  (2.12),  one  can  show  that  the  heat  flux  must  be  continuous  across  an  equi¬ 
librium  shock,  so  that  go  is  independent  of  r. 

The  system  of  governing  equations  consists  of  the  two  first-order  nonlinear 
ordinary  differential  equations  (4.4)  and  (4.5)  in  the  unknowns  u  and  6.  Because 
of  the  assumed  symmetry,  the  shock  must  be  a  cylindrical  surface  at,  say,  r  =  s. 
For  simplicity,  we  restrict  attention  to  the  case  for  which 

6a  <  06,  0«  <  0.;  (4.6) 

other  cases  can  be  treated  similarly. 


4.2  Temperature  distribution 


Under  the  restriction  (4.6)2,  Equation  (3.37)  gives 


*(7 (*0»  W)  =  M*(» r  e  («» s)-  (4-7) 

For  r  e  (a,  6),  however,  there  are  three  r mutually  exclusive  possibilities  (see  Fig. 
1  for  the  sets  Ai,  i  =  1,  2,  3) : 

(7 ,0)€Ai,  *  =  1,2,3.  (4.8) 

By  (3.37),  this  means 

*(7(r),  Hr))  =  MW),  r  €  (s,  6),  t  =  1,  2,  3.  (4.9) 

Equations(4.5)  and  (4.9)  together  with  (4.2)i  give 

r  =  b  exp[27rFj(0)/go],  r  G  [a,  6],  *  =  1,  2,  3,  (410) 

where 

Ft(B)  =  f  ki(<f>)  d<j>,  t  =  l,  2,  3.  (4.11) 

Jeh 

Since  ifc<(-),  i  =1,  2,  3,  are  positive-valued  functions,  (4.10)  can  be  inverted  for  i 
=1,  2,  3,  respectively,  to  give 

0  —  Bt(r),  re  [a, b],  i=  1,2,3.  (4.12) 

Similarly,  from  (4.5),  (4.7)  and  (4.1)i,  one  has 

r  =  o  exp{2*[Fx{6)  -  Fi(0o)]/go},  r  G  [a,  a],  (4.13) 

which  can  be  inverted  to  give 

0  =  Mr),  r  G  [a,  a].  (4.14) 

For  each  choice  of  i  in  (4.12),  continuity  of  temperature  requires 
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Oi(s )  =  60(s ),  i  =  1,  2,  3,  (4.15) 

this  leads  to  the  determination  of  the  heat  flux,  qo  =  Qi(s),  i  =  1,  2,  3. 

For  convenience,  define 


Oi(r;s,0b )  =  < 


0o(r),  re[o,s]; 


i  =  1,  2,  3.  (4.16) 

0i(r),  r  6  [s,  6], 

Here,  we  explicitly  exhibit  both  s  and  0b  as  parameters,  because  we  shall  allow 
them  to  change  in  later  chapters.  In  fact,  one  can  readily  show  that  0i(r;  s,  6b) 
is  independent  of  s. 


Note  that  in  the  case  of  0b  <  0,,  we  have  i  =  1,  and  there  is  no  possibility  of 
mixed-phase  equilibrium.  To  obtain  a  mixed-phase  equilibrium,  we  assume  that 

0b  >  0 ..  (4.17) 


Since 

fc(-,0)  =  ki(0),  far0<0 (4.18) 

we  must  have 

Oi(s;s,0b)  >0m,  t  =  2,  3.  (4.19) 

Thus  corresponding  to  each  9»,  i  =  2,  3,  s  has  a  lower  bound  r*  which  is  defined 
by 

ei(r*;r*,06)=0.,  i  =  2,  3.  (4.20) 

For  future  use,  we  define  functions  f*(-)  on  [r *  ,b]  by 

6i(f-(s);s,0b)  =  0m,  i  =  2,3.  (4.21) 


One  can  show  that 


r;(r*)  =  r;f 


rjf 


t  =  2,  3, 


(4.22) 


where  r\  is  defined  by 


(4.23) 


4.3  Displacement  field 

For,  r  €  (s,6),  there  are  three  mutually  exclusive  possibilities  corresponding 
to  (4.8),  which,  together  with  (3.30),  yield  the  strain  7(r)  as  either 


7(r)  =  r1(^,0i(r;s,06)), 

for  re  ( s,6 ), 

(4.24) 

or 

7(r)  =  r2(^,e2(r;s,06)), 

for  re  ( s,b ), 

(4.25) 

or 

7(»0  =  r3(J^,e3(r;«A)), 

for  r  e  (s,  6). 

(4.26) 

Similarly,  we  have  either 

7(r)  =ri(^:»ei(r;*,^))» 

for  re  (a,  s), 

(4.27) 

or 

7(r)  =  ri(^~»  02(r; 

for  re  (a,  s), 

(4.28) 

or 

7 (r)  =  03(r5  s,  #b)), 

/or  r  e  (a,  s). 

(4.29) 

Combining  (4.24)  with  (4.27),  we  get 

7(r)  =  ri(^;,9i(r;s,0b)), 

/or  r  G  (a,  6), 

(4.30) 

Recalling  (3.25)  and  (3.30),  we  infer  from  (4.30)  that 


-26- 


T  <  Tfi r, »,  et),  Vr  €  [rj,  f>],  (4.31) 

where 

T“(r,s,eb)  =  2irrTM(e1(r;s,0b))  =  T^(r,b}0b),  Vr  e  [rj,6].  (4.32) 

Equations  (4.30)  and  (4.2)2  give 


u(r)  =  Ui{r,a,T,0h)  =  “j^ri(^, ©!(*;«, ft))  dx,  r  e  [a, 6].  (4.33) 


If  (4.25)  and  (4.28)  hold,  we  have 


T<T2w(r,sA),  re(f5W,a], 


7T(r,  S,  eb)  <  T  <  Tef  (r,  a,  06),  r  e  [s, 6], 


«6[r5,6],  (4.34) 


where 


r2M(r,  a,  06)  =  2irrrM(G2(r;  a,  0*)), 


r  G  [r2(s),6],  s  E  [r2,6].  (4.35) 


TT(r,  a,  eb )  =  27rrrm(02(r;  a,  06)), 

Equations  (4.25)  and  (4.28),  together  with  (4.2)2,  give  the  out-of-plane  displace¬ 
ment 

u(r)  =  U2(r;s,T,0b ) 

-[  frri(^>e2(x;s,0b))dx  +  f^r2(^,e2(x;s,0b))dx],  r  E  [a,s]; 
-J?r2(£z,92(x,s}0b))dx}  r  E  [a, 6). 


{ 


(4.36) 


Finally,  the  validity  of  (4.26)  and  (4.29)  requires 

T  <  T™ (r,  s,0b),  r  E  (f|(a),a],  s  E  (r|,6]; 
T  >  T?(r,s,0b),  r  E  [a,  6],  ae[r|,6), 

where 


(4.37) 


T$* (r,  s,  0b )  =  2irrrM(G3{r\  s,  0b )) 

r  E  [f|(a),6], 

3^(r» 0d)  =  2irrrm(03(r;  a,  0*)) 


•€[r;,6].  (4.38) 
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Equations  (4.26)  and  (4.29)  together  with  (4.2)2  give 

u(r)  =  U3(r’,s,T,eb) 

-[J!ri(£z,Q3&8i0b))dx  +  J*r3(£;,e3(x;s,db))dx],  re[a,s] ; 

= 

k  -frr3(^^3(x;s,0b))dx,  r  €  [s,  6]. 

(4.39) 


4.4  Non-uniqueness  of  the  equilibrium  states 

As  we  have  seen  in  the  previous  section,  the  equilibrium  problem  under  con¬ 
sideration  does  not  have  unique  solution  for  given  T,  8a  and  Bb  in  general.  To 
examine  the  totality  of  solutions,  we  introduce  the  following  sets  in  the  s,T-plane 
for  fixed  0b: 


£i(«0  =  {(>,T)/  r<r,(«t),  «  =  &}, 

£2 (»*)  =  {(.,!■)/  8  €  [r,,6]}, 

r  T<Tf(8,eh),  a  6  (rj, (.];■) 

&(«<,)=<  (a, r)/  >. 

I  T>T?(8,eb),  a  e  [r|,  6]  J 

where 

Ti(«i)  =  mmJT1M(r,b,<>t)), 

re  [r  J  ,6] 

raM(».9t)=  irf  {^(r.a.fc)}, 

T^(s,db)=  inf  {r3M(r,s,06)}, 

re  (fi;  («),«) 

TT(a,  <W  =  max  {77>(r,  a,  *,)},  i  =  2,  3. 
re[«,b] 


(4.40) 


(4.41) 


From  (4.31)-(4.41),  we  conclude  that  a  pair  of  functions  0» (•;  s,  6b)  and 
Ui(-]s,T,9)  is  a  solution  of  our  problem  if  and  only  if  (s,T)  belongs  to  Ej(0j,),  i 
=  1,  2,  3.  Let  Qfj,  i  =  1,  2,  3,  denote  the  corresponding  three  classes  of  solutions 
as  follows 
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=  {*<(•;«, o*), Uib s,T,eb)/  (f,r)€Ei(8fc).  v*e[*.f **]}, 


(4.42) 


*  =  1,  2,  3. 

The  solutions  in  Si  axe  independent  of  the  shock  position  s;  we  will  consider  these 
solutions  as  s  =  b.  Then,  for  given  T  and  6b  appropriately,  we  have  a  solution  in 
Si,  a  one-parameter  family  of  solutions  in  S2  with  the  parameter  s,  rj  <  s  <  b, 
and  a  one-parameter  family  of  solutions  in  S3  with  the  parameter  s,  <  s  <  b. 
For  qualitative  purposes,  these  are  shown  schematically  in  Fig.  6. 


From  Fig.  6,  one  can  draw  the  following  conclusions: 


(a)  6m  <6b  <  6m: 

For  a  load  T  such  that 

0  <  T  <  T^(b,6b) 
T^ibfOb)  <T  <  Ti(8b) 

Ti(06)  <  T  <  T, 

T.  <  T  <  00 
(here  Tm  =  27rr^rM(0*)). 

(b)  6.<0b<  6*: 

For  a  load  T  such  that 

0  <  T  <  T\{9b) 


Ti(6b)  <T  <  T. 


T.  <  T  <  00 


there  exist 

a  unique  solution  which  is  smooth; 
two  one-parameter  families  of  solutions, 
one  smooth  solution; 
a  one-parameter  family  of  solutions, 
no  smooth  solutions; 
a  unique  solution; 


there  exist 

two  one-parameter  families  of  solutions, 
one  smooth  solution; 
a  one-parameter  family  of  solutions, 
no  smooth  solutions; 
a  unique  solution. 


For  sufficiently  large  load  (T  >  T„),  the  solution  involves  a  smooth  displace¬ 
ment  field  u,  but  the  temperature  field  9  is  not  smooth  unless 


*l(0.)  =  *3(0.)- 


(4.43) 


For  6*  <  6b  <0*,  the  result  is  similar  to  that  described  in  Fig.  6(a)  if  we  only 
consider  solutions  involving  at  most  a  single  shock,  or  phase  boundary.  However, 
in  the  case  6b  >  6* ,  it  may  be  necessary  to  consider  solutions  involving  multiple 
shocks. 
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5.  Quasi-static  motions  and  admissibility 
5.1  Quasi-static  motions 

Here  we  consider  quasi-static  motions  in  which,  at  each  instant,  the  state 
(0,u)  is  in  equilibrium  and  belongs  to  the  collection: 


0  =  Oi  U  U  S3-  (5.1) 

Let  L  =  [t0,  ti]  be  a  finite  time  interval,  and  let  A  be  a  mapping  of  L  into  Of;  thus 
there  are  two  functions  8(-,t)  and  u(*,<)  defined  on  [a,b]  and  such  that 

A («)  =  (*(•,*),«(•,*)),  teL.  (5.2) 

Every  mapping  A  determines  a  unique  quadruple  of  functions  *(•),«(•),  T(*)  and 
&b(-)  on  L  such  that 

*(t)  =  1,2  or  3  and  (s(t),T(t))  €  £<(t)(06(t)),  teL,  (5.3) 

and 

0(-,t)  =  9i(t)(-;s(t),0b(t)), 

t  e  L.  (5.4) 

u(-,t)  =  Vi(t)(-;s(t),T(t),0b(t)), 

The  mapping  A  is  a  quasi-static  motion  if  the  associated  quadruple  has  the 
following  properties: 

(a)  i(-)  is  piecewise  constant  on  L,  taking  one  of  three  values  1,  2,  3  at  each 
instant  t; 

(b)  s(-),  T(-)  and  9b(')  are  continuous  and  piecewise  continuously  differentiable 
on  L,  and  (s(t),T(t))  e  Ei(t)(0b(t))t  Vf  e  L; 

(c)  0(r,  •)  and  tt(r,  •)  are  both  continuous  on  L  for  every  r  e  [a,  6]. 


An  instant  t«  €  L  at  which  i(t)  is  discontinuous  is  called  a  transition  instant. 
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Proposition  5.1:  Let  *(•),  s(-),  T(*)  and  0b(-)  be  functions  on  L  with  properties 
(a)  and  (b)  above.  Define  &(-,  •)  and  tt(-,*)  on  [a,  6]  x  L  by  (5.4).  Then  6  and  u 
have  property  (c)  if  and  only  if,  at  every  transition  instant  t,  €  L, 


s(t.)  =  b.  (5.5) 

Proof:  Suppose  first  that  (c)  holds,  so  that  0(r,  •)  is  continuous  on  L.  Then  by 

(5.4),  if  t*  is  a  transition  instant  in  (to,ti), 

lim  GiW(r;s(t),0b(t))  =  lim 9i(t)(r;s(t),0b(t)),  r€[o,6].  (5.6) 

There  is  a  6  >  0  such  that  t(-)  is  continuous  on  (t*  —  tf,t.  +  6)  except  at  so 
that  (5.6)  implies 


Oj(r;s,0b)  =  Gk(r;s,0b),  r€[a,6],  (5.7) 

where  j  —  i(t*— ),  k  —  i(t.+),  s  =  s(t.),  0b  =  0b(U).  By  (4.16),  (5.7)  is 
equivalent  to 

0o(r;gj(s))  =  0o(r,qk(s)), 

0j(r’iQj(8))  =  0k(r;qk(s)), 

By  (4.5),  (5.8)  implies 

Qj(s)  =  Qk(s),  at  s  =  s(U).  (5.10) 

From  (4.5),  (5.9)  and  (5.10),  it  follows  that 

kj(0)  =  kk{0 ),  0(s(t,),t .)  <  0  <  0b.  (5.11) 


a  <  r  <  s, 

(5.8) 

8  <r  <b. 

(5.9) 

Since  t.  is  a  point  of  discontinuity  of  i(-),  one  has  j  k.  Hence  (5.11)  contradicts 
(3.38)  unless  (5.5)  holds.  Minor  modifications  of  the  above  argument  can  be  used 
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to  deal  with  the  case  t,  =  to  and  t*  =  t\.  Thus  property  (c),  in  the  presence  of 
properties  (a)  and  (b),  implies  (5.5). 

Next  we  prove  that  6(r,  •)  and  u(r,  •)  are  both  continuous  on  L  for  every 
r  €  [a,  5]  under  the  assumption  that  (5.5)  holds  at  every  transition  instant  f.  6  L. 
Note  that  continuity  of  0(r,  •)  and  u(r,  •)  at  any  non-transition  instant  •'very 
r  e  [a,  6]  follows  from  the  continuity  of  s(-),  T(-)  and  9b(-)  as  guaranteed  by 
property  (b).  At  a  transition  instant  t.  e  L,  by  property  (a)  and  (5.5),  (5.4) 
yields 


lim  (0(r,t),u(r,t))  =  (ej(r;6,0b),C^(r;6,r,06)),  re[a,b],  (5.12) 

t— 


and 


]im  (0(r1t),u(r1t))  =  (ek(r-,b,Ob),Uk(r\b,T,db)),  re  [a, 6],  (5.13) 


where  j  =  i(i*— ),  k  =  t(t,+)  and  T  =  T(t*),  6b  —  6b(t *).  By  the  analysis  in  the 
previous  chapter,  one  has 


(©*(r;  b,  &b),  Ui(r;  6, T,  0&))  =  (©1(r;6,^),(/i(r;6,r,tfb)) 

(5.14) 

re[a,b],  i  =  l,  2,  3. 

Hence  0(r,  •)  and  it(r,  •)  are  continuous  at  a  transition  instant  and  therefore  at 
all  instants  in  L.  Thus  (5.5),  in  the  presence  of  properties  (a)  and  (b),  implies 
property  (c),  and  the  proposition  is  established. 


5.2  Admissibility 

For  quasi-static  motions  in  the  annular  cylinder  problem,  Equation  (2.20) 
which  is  a  consequence  of  the  second  ’aw  of  thermodynamics  reduces  to1 


1  Here  the  normal  N  is  chosen  towards  the  inner  wall  of  the  annular  cylinder. 
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-/(s(t),t)s(i)  >0,  t  e  L,  (5.15) 

where  the  driving  traction  f  is  given  by 

/(s,t)  =fci(t)(^,0(M)),  (5.16) 

and  hj(-,  •)  are  the  material  functions  defined  by 


hj(r,  9)  ='J,(r1(r,  0),  9)  -  ^(t,  9),  9)  -  r^r,  9)  -  T,(t ,  0)], 
/or  rm(0)  <  t  <  ta/(5),  0.  <  0  <  0*,  j  -  1,  2,  3, 

where  is  defined  by  (3.21). 


(5.17) 


We  will  say  a  quasi-static  motion  is  admissible  if  (5.15)  holds. 


By  the  definitions  (5.17)  and  (3.30)  of  hj  and  1^,  one  can  see  that  hi(-,  ) 
vanishes  identically,  and  that 


dhjjr,  9) 
dr 


j  —  2,  3, 


(5.18) 


for  Tm(9)  <  t  <  tm(0),  9  e  (0.,  9*). 

Recalling  the  definition  in  Chapter  3  of  the  Maxwell  stress  rc(0)  and  noting  that 
^(^Af (0), 0)  =  ri(rM(0), 0)  =  7Af(0)  for  0  €  [0*,  0*],  one  can  show  that 


/•tm  (®) 

Mr,  0)  =  -  J  [r2(o,  0)  -  rx((7, 0)]  do, 

Mr,  0)  =  /r  [r3(<r,  0)  -  I^o,  0)]  do, 

Jr c{9) 

for  Tm(0)  <r  <tm(9),  9  6  (0*, 0*). 


(5.19) 


Therefore,  we  have  the  following  proposition: 

Proposition  5.2:  Let  A  be  a  quasi-static  motion  on  L.  If,  at  some  instant 

t  6  L,  A(<)  €  $>2,  then  the  admissibility  condition  (5.15)  requires  that  either 
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i(t)  >  0, 

(5.20) 

/(«{(),  t)  =  0. 

(5.21) 

This  proposition  follows  immediately  from  (5.15),  (5.16)  and  (5.19)i. 

As  a  consequence  of  the  admissibility,  the  next  proposition  excludes  the  pos¬ 
sibility  that  a  particle  transforms  to  phase  2  from  either  phase  1  or  phase  3. 

Proposition  5.3:  Let  A  be  an  admissible  quasi-static  motion  on  L.  If  at  some 
instant  £'  €  [£0,£i),  A(£')  e  QfiUS>3,  then  A(£)  €  9iU93  for  any  instant  £  €  [£',£i]. 

Proof:  Suppose  the  result  of  the  proposition  is  false.  Then  there  exists  an 

admissible  quasi-static  motion  A  on  L  such  that 


A(t')  €  3i  U  33,  A(t")  e  9-2,  (5.22) 

where 

£'  <£",  [t',t”]cL.  (5.23) 

Since  A (£")  €  $2  implies  that  s(£")  <  6,  by  Proposition  5.1,  there  exists  a  finite 
interval  [£*,<2]  C  [£',  £"]  where  £,  is  a  transition  instant  such  that 

r  i(£)  <  0,  for  £.  <  £  <  £2, 

<  (5.24) 

A(£)  e  9?2,  for  £*  <  £  <  £2. 

Hence  Proposition  5.2  requires  that  the  driving  traction  /(«(•)»  •)  vanish  identi¬ 
cally  on  [£.,£2]  which,  by  (5.16)  and  (5.19)i,  implies  that 

T(t) 

2tt s(t)  =  ta*(^(s(^)i  0)»  Vt  ^  [£*,£2], 
i.e. 

ri(^)’eW!),t))  =  rs(^|)’#{,(t)’‘))  = 

for  t  e  [£.,£2]. 


(5.25) 

(5.26) 
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Noting  that 


gfCr. 9) 

dl 


=  0,  far  7  =  7 m(0),  0  €  [0*,  5*], 


(5.27) 


one  finds  from  (4.4)  and  (4.5)  that 

dd  dd( 

^>-aF(r’‘)“^f)+a?(r>*)' 

lim  [/l(»(r,t))^(r,«)]=  lim j/j(«(r,t))j{r,t}|, 

r— C7T  r— >*(t)+  OT 

far  t*  <t  <  t2- 

By  (3.38),  Equation(5.28)  holds  only  if 


(5.28) 


which  implies  that 


0(s(t),t)  =  0.,  vte[t.,t2], 


(5.29) 


0.  =  lim  0(s(t),t)  =  0(6,t,)  =  0b(t.).  (5.30) 

This  contradicts  the  assumption  made  in  (4.17),  establishing  the  proposition. 


According  to  Proposition  5.3,  an  admissible  quasi-static  motion  that  begins 
in  either  phase  one  ($i)  or  phase  three  ($3)  can  never  enter  phase  two  ($2). 
In  particular,  quasi-static  motions  that  commence  at  the  undefonned  state  can 
never  enter  phase  two.  We  thus  exclude  states  in  $2  from  here  on.  Recalling  the 
definitions  (4.16),  (4.33)  and  (4.39)  of  0*,  U\  and  U 3,  one  can  show  that 


[0i(-;6,^),^i(-;^r,^)]  =  [e3(-;6,^),^(-;6,r,<?6)], 

far  T  <  Ti(9b). 

Noting  the  definitions  (4.40)  and  (4.42)  of  fi(6b)  and  one  can  write 


(5.31) 


Ol  =  {«3  h)M-,»,T,9t)/  (..TJtLift)),  (5.32) 

Thus,  for  any  quasi-static  motion  that  begins  in  either  phase  one  (Si)  or  phase 
three  (§3),  we  have 
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6(r,t)  =  e3(r;s(t),0b(t)), 

V(r,t)  e  [a,b]  x  L.  (5.33) 

u(r,  t)  =  U3(r\ s(t ),  T(t),  0b(t)), 

Hence,  when  the  load  T(t)  and  the  temperature  0b(t)  on  the  outer  wall  are 
given,  the  temperature  field  6(-,  t)  and  the  displacement  field  u(-,t)  for  each  in¬ 
stant  t  €  L  are  determined  by  (5.33)  uniquely  up  to  an  unknown  function  s(-) 
which  specifies  the  location  of  the  shock,  or  phase  boundary,  at  every  instant.  As 
pointed  out  by  Abeyaratne  and  Knowles  [13],  this  suggests  that  some  additional 
constitutive  information  concerning  the  kinetics  of  the  phase  transformation  must 
be  provided. 

For  the  so-called  hard  device  problem,  the  displacement  £(•)  on  the  inner 
wall  would  be  specified  instead  of  the  load  T(-).  By  (5.33), 

6  =  6(T-,  s,  0b)  =  U3(a\  s,  T,  0b),  at  any  instant.  (5.34) 

Recalling  (4.39),  (3.30)  and  (3.22),  one  can  show  that 

A 

J^(7V,  )<0,  tor  all  T,  (5.35) 

assuring  that  there  is  a  unique  inverse  T(-’,s,0b)  for  6(-;s,0b ): 


T  =  f(6;s,0b) 

for  all  6 ,  and  for  appropriate  values  of  s  and  Ob- 
Define 

U(r;s,6,0b)  ~  U3(r-,s,f(S\s,0b),0b), 

for  r  e  [a,  6],  s  6  [r|,  6],  6  >  0,  0b  >  0a • 

Then  by  (5.33), 

0(r,t)  =  93(r;s(t),0b(t)), 
u(r,t)  =  U(r;s(t),6(t),0b(t)), 


(5.36) 


(5.37) 


V  (r,  t )  £  [a,  b]  x  L. 


(5.38) 


-37- 


Again,  the  quasi-static  motion  is  described  by  (5.38)  uniquely  up  to  an  unknown 
function  s(-). 
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6.  Completion  of  constitutive  assumptions 

As  we  have  seen  in  Chapter  5,  additional  constitutive  information  is  needed 
in  order  to  determine  the  location  of  the  phase  boundary  during  quasi-static  mo¬ 
tions.  One  may  view  this  constitutive  deficiency  as  arising  from  the  fact  that 
the  constitutive  relations  (2.21)  -  (2.24)  are  not  meaningful  on  the  phase  bound¬ 
ary.  In  this  chapter,  we  complete  the  constitutive  description  of  the  material  by 
specifying  a  kinetic  relation  that  controls  the  rate  at  which  the  phase  transfor¬ 
mation  takes  place,  and  a  transformation  strain  that  determines  the  onset  of  the 
transformation. 

6.1  The  kinetic  relation 

As  in  the  work  of  Abeyaratne  and  Knowles  [13]  on  modeling  phase  transfor¬ 
mations  in  the  setting  of  finite  elasticity,  we  assume  that  the  normal  component 
of  the  propagation  velocity  of  phase  boundary  depend  on  the  driving  traction 
and  the  temperature  on  the  phase  boundary.  Thus  we  suppose  that  there  exists 
a  function  V{-,  6)  for  6  €  such  that 

V)v(x,<)  =  y(/(x,  t),0(x,  t)),  for  x  6  S(t),  t  e  L.  (6.1) 

For  the  annular  cylinder  problem,  we  have  shown,  in  the  previous  chapter, 

that 

/(s,t)  =  for  (s,t)  €  [a, 6]  x  L,  (6.2) 

where  h( •,  •)  is  a  material  function  defined  through  the  inverse  functions  T3,  r i 
of  (3.30)  by 

h(r ,0)=r  [r3(<r,0)  -  r,(<r,0)]  cfcr. 

Jr,(S) 

forTm(0)  <t<tm(0),  fle[0.,0']t 


(6.3) 
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where  rc{6)  is  the  Maxwell  stress  at  temperature  9. 

Define 

/»(*)  =  Krm(9),9)t  fxi(0)  =  Ktm(9),9),  far  9  G  [9.X]-  (6.4) 

Then  by  (6.3), 

M9)  <  h(r}9)  <  MB ),  for  rm(9)  <  r  <  tm(9),  9  G  [*.,*•].  (6.5) 

Hence  in  the  cylinder  problem,  the  kinetic  response  function  y(-,  9)  is  to  be  defined 
on  (/m(0)»/Af(0))  for  9  G  [0*,0*],  and  the  kinetic  relation  (6.1)  specializes  to 

-s  =  V (/(s,  t),  0(s,  t)),  for  (s,  t )  G  [a,  6]  x  L,  (6.6) 

where  s(t)  is  the  radius  of  the  phase  boundary  at  time  t. 

Concerning  the  form  of  the  kinetic  function  (see  Fig.7),  we  assume  that  V (•,  9) 
is  continuous  on  (fm(9),  fxf(9))  for  every  0  G  (0*,0*)  and  that 

(i)  There  exist  two  material  functions  m{9)  and  M{9)  on  (9»,9*)  satisfying 

fm(9)  <  m{9)  <  0  <  M{9)  <  fM(B),  far  9  G  (0.,  **),  (6.7) 

and  such  that 

1<  0  for  fm(9)  <  f  <  m(9); 

=  0  for  m(9)  <f<  M[9)\  9  G  ( 9., 9 *),  (6.8) 

>  0  for  M{9)  <f<  f„(9); 

(ii)  moreover 


V(O,0)  =  O, 


for  9  G  [0*,0*], 


(6.9) 
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lim  V(f,  0)  =  — oo,  lim  V(/,0)  =  oo,  /or  0  G  (0.,0*).  (6.10) 

Note  that  m(0)  =  M{6)  =  0  is  a  permissible  choice. 

If  T(t)  is  given  and  6(s{t),t)  is  represented  in  terms  of  s(t)  by  (5.38)  and 
(4.16),  then  (6.6)  become  an  ordinary  differential  equation  for  s(t).  To  arrive  at 
an  appropriate  initial  condition  for  this  differential  equation,  we  need  a  criterion 
for  initiation  of  the  phase  transformation,  to  which  we  now  turn. 


6.2  The  transformation  strain  postulate 

Define  two  functions  7C(*)  and  7C(-)  on  [0*,  0*]  by 


1.(0)  =  IMt.P),  0),  y=(0)  =  r3(rc(0),  0), 
where  tc(0)  is  the  Maxwell  stress. 


(6.11) 


Assume  that  there  exist  two  material  functions  yr(m)  and  7T(-)  on  (0,,  0*), 
taking  values  (see  Fig.  8) 


7c (0)  <  7r(0)  <  7Af(0)  <  7m(0)  <  7T(0)  <  7C(0)> 

for  0  E  (0*,  0*), 


(6.12) 


such  that 


i)  A  particle  at  r  =  f  will  “spontaneously”  transform  from  phase-1  to  phase-3 
at  t  =  t  if 


0(f,t)  G  (0.,0*)  and  7(r,f)  =  7r(0(r,t));  (6.13) 

ii)  A  particle  at  r  =  f  will  “spontaneously”  transform  from  phase-3  to  phase- 1 
at  t  =  t  if 
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0(r,t)  €  (0.,00)  and  n/(f,i)  =  'yT(6(f,i));  (6.14) 


We  may  call  7t(0  and  7T(')  the  (1,3)- transformation  strain  and  the  (3,1)- 
transformation  strain,  respectively. 
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7.  Macroscopic  thermo- mechanical  responses 

We  now  elucidate  the  features  of  the  macroscopic  response  of  the  cylinder 
under  various  loading  programs,  with  the  temperature  6b  at  the  outer  wall  fixed. 
We  shall  adopt  a  kinetic  relation  of  the  type  illustrated  in  Fig.  7.  For  convenience, 
we  recall  here  from  (5.34)  and  (5.36)  the  macroscopic  response  for  soft  and  hard 
devices,  respectively: 


6(t)  =  6(T(t);  a(t),  db(t)),  t0<t<  tl5  (7.1) 

T(t )  =  f  (*(i);  s(t),  6b(t)),  t0<t<  ti.  (7.2) 


7.1  Hysteresis  and  rate-dependence 

First  let  T(t)  —  At,  corresponding  to  soft  device  loading  at  a  constant  rate 
A  from  the  undeformed  state.  Assume  that  the  final  value  of  the  load  T  is  larger 
than  T.(0b),  and  that  the  fixed  value  of  the  outer  wall  temperature  6b  is  such 
that  6+  <  6b  <  0«.  The  T  —  6  response  is  shown  schematically  in  Fig.  9.1  After 
loading  begins,  the  point  rises  from  the  origin  O  along  the  response 

curve  OAQ,  associated  with  a  single  phase  solution  at  each  instant.  At  t  =  t, 
the  point  ( 6(t),T(t ))  reaches  the  point  Q  in  Fig.  9.  When  this  occurs,  the 
amount  of  shear  7(6,  t)  at  the  outer  wall  achieves  the  value  'yr(6b)  of  the  (1,3)- 
transformation  strain  at  the  wall  temperature  6b,  so  that  a  phase  transformation 
is  initiated  according  to  the  initiation  criterion  (6.13),  and  a  1,3-phase  boundary 
emerges  at  r  =  b.  The  kinetic  relation  (6.6)  then  takes  over,  and  the  phase 
boundary  has  the  initial  velocity  s(t)  =  —V(f(b)  t),  6b).  Under  the  control  of 
(6.6),  8(t)  monotonically  decreases  with  t  until  s  —  r£(6b);  this  occurs  at  time  t  = 
T*(6b)/ A.  The  axial  displacement  6(t)  at  the  inner  wall  is  given  by  (7.1),  and  point 

1  Here  it  is  assumed  that  0*  <  6b  <  this  guarantees  the  absence  of  residual 
deformation,  which  is  considered  in  the  next  section. 
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(6(t),T(t))  moves  along  the  curve  QBP.  This  stage  of  the  process  is  accompanied 
by  a  phase  transformation.  As  the  load  T(t)  continues  to  increase,  the  phase 
boundary  remains  at  r  =  rj(0j,)  since  the  temperatures  at  points  for  which  r  < 
73 (6b)  are  too  low  (  <  0.)  for  the  material  to  support  a  phase  transformation. 

Now  suppose  that  T(t)  is  decreased  at  the  constant  rate  A  from  its  largest 
value  at  the  point  Z  in  Fig.  9  to  zero.  The  associated  response  curve  at  first 
follows  the  arc  ZP,  and  the  phase  boundary  stays  at  r  =  (#{,),  so  there  is  no 

reverse  phase  transformation  at  this  stage.  As  T(t)  continues  to  decrease,  s(t) 
increases  with  time  under  the  control  of  (6.6),  i.e.,  the  phase  boundary  moves 
toward  the  outer  wall,  and  the  point  ( S(t ),  T(t))  moves  along  the  curve  PCA.  As 
T(t)  decreases  to  zero  from  its  value  at  A,  s(t)  remains  at  the  value  b  until  the 
point  (6(t),T(t))  returns  to  the  origin  along  the  curve  AO. 

If  we  change  the  loading  rate  A,  the  arcs  QBP,  PCA  will  change,  but  OAQ 
and  PZ  would  remain  the  same.  For  different  loading-unloading-reloading  pro¬ 
cesses,  we  have  different  macroscopic  response  curves  similar  to  those  predicted 
by  Abeyaratne  and  Knowles[13]. 

7.2  Residual  defomation 

If  the  outer  wall  temperature  0&  is  fixed  such  that  0,  <  Of,  <  6*,  some  residual 
deformation  would  be  produced  by  the  mechanical  cycling  process  described  in 
the  previous  section.  In  this  case,  the  macroscopic  response,  during  the  loading 
period,  is  qualitatively  the  same  as  that  sketched  in  the  last  section,  but  during 
the  unloading  period,  it  is  quite  different.  At  first  the  response  curve  (see  Fig. 
10)  follows  the  arc  ZP,  and  the  phase  boundary  remains  stationary  at  r  =  7-3  (0j,). 
When  T(t)  decreases  to  zero  from  T,(0&),  according  to  (6.6),  the  phase  boundary 
moves  towards  the  outer  wall  and  arrives  at  r  =  si  <  6,  corresponding  to  the  arc 
PCD.  If  the  driving  traction  f  at  the  phase  boundary  at  this  instant,  say  t  =  ti, 
is  less  than  the  critical  value  m(0(si,fi)),  then  the  phase  boundary  continues  to 
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move  towards  the  outer  wall,  and  the  point  (6(t),  T(t))  “creeps”  to  E  from  D  at 
zero  load.  Creep  response  under  more  general  conditions  will  be  discussed  in  the 
next  chapter. 

7.3  Stress  relaxation 

Let  the  outer  wall  temperature  Ob  be  fixed,  with  0»  <  Ob  <  0* .  We  first  let 
6(t)  =  At,  corresponding  to  hard  device  loading  at  a  constant  rate  A  from  the 
undeformed  state.  The  corresponding  response  curve1  is  the  arc  OQB  in  Fig.  11. 
If  we  stop  loading  at  B,  subsequently  keeping  6  constant,  the  phase  boundary 
continues  to  move  towards  the  inner  wall,  and  the  point  (S(t),  T(t ))  approaches 
the  point  C  at  which  the  driving  traction  is  equal  to  the  critical  value. 


It  is  assumed  in  Fig.  11  that  9*  <0b  <  0,. 


l 
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8.  Creep  due  to  phase  transformation 


As  described  in  Chapter  4,  the  experiment  of  Sammis  and  Dein[14]  suggests 
that  a  solid  may  creep  at  constant  load  when  it  experiences  a  phase  transforma¬ 
tion.  In  this  chapter,  we  discuss  the  possible  occurrence  of  creep  as  predicted  by 
our  model,  and  we  show  that  the  predictions  are  in  qualitative  agreement  with 
the  experimental  results  reported  in  [14]. 


We  consider  loading  by  a  soft  device,  with  the  load  T(t)  at  t  given  by 


f  (t)  =  l 


j'  Tt/tu  0  <  t  <  ti; 
[  T,  t  >  t,. 


(8.1) 


corresponding  to  ramp  loading  until  time  ti,  followed  by  a  constant  load  T.  The 
outer  wall  temperature  6b  is  kept  constant  in  the  range  (0*,  0*).  The  annular 
cylinder  is  said  to  creep  if  the  displacement  6{i)  at  the  inner  wall  varies  with  time 
in  the  time  interval  [<i,  oc)  during  which  the  load  remains  constant  at  the  value 
T.  After  the  constant-load  stage  is  reached,  the  creep  rate  at  the  inner  wall  is 
found  from  (7.1)  to  be  given  by 


•  KJ  O 

6(t)  =  —  (T;a(t),eb)8{t),  for  t  >  tu  (8.2) 

where  the  radius  s(t)  of  the  phase  boundary  at  time  t  is  the  solution  of  the 
ordinary  differential  equation  (6.6)  that  satisfies  the  initial  condition 


s(h)  =  si.  (8.3) 

Here  Si  €  (r^,  6]  is  the  radius  of  the  phase  boundary  at  the  end  of  the  ramp 
loading  period.  One  may  have  either  si  =  b  or  si  <  b.  In  the  former  case,  the 
ultimate  load  level  T  is  too  small  to  have  initiated  a  phase  transformation  during 
the  ramp  period.  In  the  latter  case,  a  phase  transformation  has  been  initiated 
during  the  ramp  period,  and  the  phase  boundary  has  already  moved  inward  to 
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r  =  si  by  the  time  the  ultimate  load  level  T  is  reached.  (See  the  discussion  in 
the  previous  chapter  concerning  constant-rate  loading  processes.) 

Since  by  (5.33)  the  temperature  distribution  9{r ,  t )  depends  on  the  time  only 
through  the  phase  boundary  radius  s(t)  when  the  outer  wall  temperature  6b  is 
kept  constant,  one  may  write  6{s{t),t)  =  0(s(t)f0b),  where  we  have  also  indicated 
explicitly  in  the  notation  the  dependence  on  the  outer  wall  temperature  9b  ■  By 
(6.2),  we  may  then  rewrite  the  kinetic  relation  (6.6)  in  the  form 

-s(t)  =  V(s(t),T,0b),  for  t  >  ti,  (8.4) 

where 

V(s-,T,eh)  =  v(h(~,6(s,eb)),e(s,e »»,  (8.5) 

while  V  is  defined  by  (6.1).  We  note  that  the  differential  equation  (8.4)  is  au¬ 
tonomous. 

According  to  (8.2),  the  creep  rate  6(t)  is  zero  if  the  radius  of  the  phase 
boundary  s(t)  is  constant  for  t  >  t\.  This  occurs  if  and  only  if  the  differential 
equation  (8.4)  has  the  solution 


s(t)  =  si,  for  t  >  t\. 
This  in  turn  happens  if  and  only  if 


(8.6) 


V(»i;T,ft)-0.  (8.7) 

Thus  creep  is  absent  if  and  only  if  the  initial  radius  si  is  a  zero  of  the  kinetic 
response  function  V*(-]T,0b).  Whether  or  not  (8.7)  holds  will  depend  not  only 
on  the  ultimate  load  level  T  and  the  value  of  the  outer  wall  temperature  9b ,  but 
in  general  on  the  particular  kinetic  response  function  and  on  the  detailed  values 
of  the  other  constitutive  parameters  of  the  model  as  well. 


-47- 


For  the  loading  program  (8.1),  for  a  specific  kinetic  response  function  V  of 
the  type  shown  in  Fig.  7,  and  for  suitably  chosen  values  of  the  various  constitutive 
parameters,  we  have  numerically  solved  the  initial  value  problem  (8.4),  (8.3)  and 
determined  the  creep  rate  according  to  (8.2)  for  four  different  values  0i,  62,  #3,  #4 
of  the  outer  wall  temperature,  ordered  by  <  62  <  O3  <  #4.  The  results  are 
shown  in  Fig.  12  as  plots  of  6  vs.  time  t.  The  computation  shows  that  no  creep 
occurs  at  the  two  lower  temperatures,  but  significant  creep  rates  have  appeared 
at  the  higher  temperatures.  The  qualitative  features  of  these  results  agree  with 
the  observations  of  Sammis  and  Dein  [14],  which  are  sketched  schematically  in 
Fig.  13. 

Our  numerical  results  in  Fig.  12  seem  to  suggest  that 


lim  6(t)  =  0.  (8.8) 

t— *00 

To  see  this  long-term  feature  of  creep,  we  assume  that  creep  occurs  after  the  ramp 
loading  period  under  consideration.  Let  the  arc  OQA  in  Fig.  14  represent  the 
response  curve  corresponding  to  the  ramp  loading  period  [0,ti].  The  subsequent 
response  must  correspond  to  a  horizontal  line  segment  starting  at  the  point  A, 
say  AB.  By  the  property  of  the  autonomous  equation  (8.4),  the  final  radius  sj  of 
the  phase  boundary  corresponding  to  the  point  B  must  be  such  that 


v(«;r, 


>0, 
=  0, 


s  >  s 

S  =  S  /. 


This  assures  that  the  autonomous  equation  (8.4)  can  be  integrated  as 


(8.9) 


tn  dr 

-w.  W5'  s>s/- 

The  integral  in  (8.10)  is  divergent  as  s  tends  to  s/  if 


(8.10) 


lim 

»-»/ 


V(»;r,gt) 

8  —  8} 


=  0. 


(8.11) 
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Any  kinetic  response  function  V*  with  the  property  (8.11)  thus  gives 


lim  s(t)  =  «/,  and  lim  s(t)  =  0.  (8.12) 

t— ♦  oo  t-^oo 

The  long-time  result  (8.8)  then  follows  from  (8.2)  and  (8.12). 
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9.  The  shape-memory  effect 

For  9  e  (#.,<?*),  the  Helmholtz  free  energy  potential  \&(-,  9)  in  (3.21)  has  two 
wells  (see  Fig.  15),  with  corresponding  local  minima  at  7  =  0  and  at  7  =  70 (9)  and 
correspondingly,  the  stress-strain  response  curve  has  the  form  shown  in  the  Fig. 
3(c).  This  allows  the  body  to  achieve  an  unstressed  but  deformed  state  through  a 
mechanical  cycling  process  at  the  fixed  temperature  9 ,  starting  at  the  undeformed 
state.  The  residual  deformation  can  be  removed  through  a  temperature  cycling 
process,  because  the  unstressed  state  must  also  be  undeformed  as  long  a s  0  € 
(0,0.)  U  (0*,  00).  Fig.  16  illustrates  the  so-called  shape-memory  behavior. 


To  see  the  shape-memory  effect  in  the  simple  setting  of  finite  anti-plane  shear 
of  the  annular  cylinder,  we  consider  a  mechanical  cycling  process  T  =  T(t),  6b  = 
constant,  followed  by  a  thermal  cycling  process  9b  =  9b(t),  T  =  constant: 


’T0t,  0  <  t  <  1; 

T(t)  =  {  To (2  —  t),  1  <t<2; 

10,  2  <  t  <  3, 


eb  =  eb(t)  =  { 


0i  +  f  (#1  —  ^*)(2  —  t), 
l*! -5(0i-0.)(3  -0, 

where  6\  €  (0*,0*)  and  To  >  T*(0i). 


0i » 


0  <  t  <  2; 

2  <  t  <  2.8; 
2.8  <  t  <  3, 


(9.1) 


(9.2) 


The  macroscopic  response  arising  from  this  program  of  thermo-mechanical 
loading  is  shown  schematically  in  Fig.  17.  During  the  mechanical  cycling  period, 
the  T  —  6  response  is  represented  by  the  path  OQPZD,  and  some  residual  defor¬ 
mation  is  clearly  produced  (see  Section  7.2  for  details).  The  mechanical  load  T 
vanishes  identically  during  the  temperature  cycling  period;  the  thermal  response 
6  vs.  6b  is  shown  in  Fig.  17(b).  After  the  wall  temperature  6b  begins  to  drop,  the 
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phase  boundary  moves  towards  the  outer  wall,  and  the  point  (6(t),9b(t))  moves 
along  the  curve  DE.  At  the  point  E,  the  phase  boundary  disappeared  at  the  outer 
wall  and  the  residual  deformation  is  removed.  The  point  (6(t),9b(t))  moves  along 
the  line  segment  EF  as  6b  continues  decreasing  to  9m ,  and  it  moves  along  the  line 
segment  FO  when  6b  increases  from  9,  to  9\.  Neither  phase  transformation  nor 
deformation  accompanies  this  stage  of  the  process. 

We  have  made  a  numerical  simulation  to  demonstrate  the  shape-memory 
behavior  of  the  annular  cylinder  during  the  thermo-mechanical  loading  process 
described  by  (9.1)  and  (9.2).  For  appropriately  chosen  values  of  the  constitutive 
parameters,  Fig.  18  shows  the  configurations  of  half  a  cross  section  of  the  body 
at  a  sequence  of  time  instants. 

During  the  mechanical  cycling  process,  a  phase  transformation  is  initiated 
and  a  1,3-phase  boundary  emerges  at  the  outer  wall  of  the  cylinder  when  the 
mechanical  load  reaches  a  certain  level.  The  phase  boundary  moves  towards  the 
inner  wall  until  it  coincides  with  the  surface  r  =  r^fli)  as  the  load  continues  to 
increase.  This  is  shown  by  Fig.  18(a)  and  (b)  in  which  the  darker  area  represents 
the  part  of  the  body  where  the  material  is  in  phase-3.  As  we  have  discussed 
in  Chapter  7,  the  phase  boundary  remains  at  r  =  r£(9i)  when  the  load  T(t) 
continues  to  increase  (see  Fig.  18(b)  and  (c)).  From  Fig.  18(c),  (d)  and  (e),  one 
can  see  that  the  load  T(t)  decreases  from  its  maximum  To,  the  phase  boundary 
stays  at  r  =  rl(9i),  and  it  moves  backwards  only  after  the  load  drops  below  a 
certain  level.  Fig.  18(f)  shows  that  at  the  end  of  mechanical  cycling  process, 
the  body  is  in  a  two-phase  equilibrium  state  and  some  residual  deformation  has 
been  produced.  During  the  temperature  cycling  process,  the  phase  boundary 
continues  to  move  towards  the  outer  wall  and  disappears  there  as  the  outer  wall 
temperature  decreases.  This  is  shown  by  Fig.  18(f),  (g)  and  (h).  Finally,  Fig. 
18(h)  and  (i)  show  that  the  body  remains  at  the  undeformed-unstressed  state 
during  the  stage  of  temperature  increasing.  Therefore,  we  see  that  the  residual 
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deformation  produced  by  the  mechanical  cycling  process  has  been  removed  during 
the  temperature  cycling  process  and  the  body  has  achived  its  original  shape. 
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Appendix.*  A  note  on  the  constitutive  parameters 


As  mentioned  in  Chapters  8  and  9,  we  have  done  some  numerical  work  to 
exhibit  the  creep  feature  and  the  shape-memory  behavior  predicted  by  our  model. 
It  should  be  mentioned  that  the  available  experimental  data  is  insufficient  to  permit 
the  determination  of  all  constitutive  parameters  needed  in  the  model  constructed 
here.  For  the  purpose  of  making  qualitative  comparison  with  experiments  we  have 
made  rather  arbitrary  choices  of  these  parameters.  For  those  readers  who  may  be 
interested,  we  give  here  some  details  regarding  our  computation. 


We  have  assumed  that  the  heat  conductives  k t  and  fc3  are  two  different  con¬ 
stants,  and  we  have  approximated  the  kinetic  response  function  V  by  a  piecewise 
linear  function: 


V(f,  8)  =  | 


»m  (/  - 

o, 

M/  ~ 


&m/m  <  f  <  fur, 

<f<  V0  6  [0.,  (/*). 

fm(0)  <  /  <  rfm/m(0), 


where  s*f,sm,dAf  and  dm  are  all  constants. 


The  numbers  we  assigned  tr  the  constitutive  parameters  are  listed  in  the  fol¬ 
lowing  table;  the  first  row  was  used  for  modeling  creep  and  the  second  was  used  for 
modeling  shape-memory  behavior. 


k3/ki 

a 

00  /« 

n 

m 

s,„ 

SM 

dm 

mm 

1.5 

.04 

1.35 

3 

9.5 

.001 

.01 

2  x  10fi 

4  x  107 

1.5 

.04 

1.23 

3 

2.3 

.001 

.3 

2  x  10* 

3  x  10* 
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One  may  notice  that  the  values  of  sm  and  d, m  in  the  first  row  are  much  smaller 
than  those  in  the  second  row;  this  results  in  the  fact  that  the  phase  boundary 
propagates  comparatively  more  slowly,  so  that  we  have  a  better  view  in  the  plots 
of  creep  rate  vs.  time  shown  in  Fig.  12. 
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Fig.  3.  Stress  strain  response  curves 
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Fig.  7.  Kinetic  response  function 


Fig.  9.  Loading-unloading  path 
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Fig.  10.  Residual  deformation 


Fig.  11.  Load  relaxation 
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Fig.  12.  Plots  of  creep  rate  8  vs.  time  t  as  predicted  by  the  present  model 
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Fig.  17-  Macroscopic  response  of  the  cylinder  during  a  process 
e  illustrating  shape-memory  effect 


Fig.  18.  Selected  frames  from  simulation  of  shape-memory  effect 
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Fig.  18.  Selected  frames  from  simulation  of  shape-memory  effect 
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